The freeze-out mechanism and phase-space density in 
ultrarelativistic heavy-ion colhsions 
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We explore the consequences of a freeze-out criterion for heavy-ion collisions, based on pion escape 
probabilities from the hot and dense but rapidly expanding collision region. The influence of the 
expansion and the scattering rate on the escape probability is studied. The temperature dependence 
of this scattering rate favors a low freeze-out temperature of ~100 MeV. In general, our results 
support freeze-out along finite four-volumes rather than sharp three-dimensional hypersurfaces, 
with high-p_L particles decoupling earlier from smaller volumes. We compare our approach to the 
proposed universal freeze-out criteria using the pion phase-space density and its mean free path. 



I. INTRODUCTION 

Experiments with nuclear collisions at ultrarelativistic 
energies are performed with the aim of studying strongly 
interacting matter under high temperature and density. 
Such conditions are created very early in the collision. 
Hadrons at low transverse momentum form the bulk of 
the particles produced in these collisions. They are pro- 
duced at the end of the hot and dense partonic phase, but 
they subsequently scatter in the confined hadronic phase 
prior to decoupling from the collision system ("freeze- 
out" ) . In general, the collective evolution of the hot and 
dense matter leaves a distinct imprint on the phase space 
distribution of the fireball at freeze-out. For example, 
pressure gradients generated in the early stage of the col- 
lision and indicative of the equation of state result in col- 
lective transverse, radial and elliptic flow which shape im- 
portant features of hadronic one- and two-particle spec- 
tra. However, to disentangle such information from fea- 
tures generated during freeze-out a refined understanding 
of the decoupling process is needed. This is the condi- 
tion which has to be satisfied locally for decoupling to 
take place. 

Freeze-out is often modeled as a sudden breakup of 
the fireball on a fixed three-dimensional hypersurface , 
specified e.g. by a critical freeze-out temperature or den- 
sity. This implements in a very simplified way the phys- 
ical picture that the hadronic scattering rate drops with 
particle phase space density and determines freeze-out. 
However, the current difficulties in interpreting the HBT 
spectra measured at RHIC Ij motivate to go beyond 
a sharp three-dimensional implementation of freeze-out. 
Moreover, the recently observed increase in average pion 
phase space density from SPS {^/s = 17 AGeV) to RHIC 
(-y/s = ISOAGeV) (SltlQi indicates that — in contrast to 
an earlier suggestion — additional factors beyond the 
particle phase space density must be taken into account. 
This prompts us in the present work to discuss quantita- 
tively how hadrochemical composition [aH, local tem- 
perature and collective velocity gradients '5'| , and particle 
momentum influence the decoupling of particles from the 
collision region. Since hadrochemistry does, and freeze- 
out temperature or velocity gradients may change in go- 



ing from SPS to RHIC energies, this study also addresses 
the possible dependence of hadronic freeze-out on bom- 
barding energy. 

Clearly, in the general case, particles decouple from the 
fireball gradually [^Q, as it is naturally implemented e.g. 
in cascade generator simulations. The question arises to 
what extent this generic behavior can be described by 
the limiting case of a sharp three-dimensional freeze-out. 
To discuss this question, we first introduce in Section 
IhI the so-called escape probability which characterizes 
the probability of a particle to decouple from the sys- 
tem. In Section IlIII we demonstrate on simple examples 
how the escape probability is related to the expansion 
strength and the scattering rate. This illustrates the old 
suggestion |^ , that freeze-out happens when the dilution 
rate becomes comparable or larger than the scattering 
rate. Moreover, the scattering rate depends on densi- 
ties of individual species and the corresponding cross- 
sections 0,|3|- For pions, e.g., scattering on nucleons is 
much more important than scattering on other pions due 
to the larger cross-section of the former. In Section HVI 
we discuss the impact of the hadrochemical composition 
on the scattering rates and freeze-out at SPS and RHIC. 
This also allows to address the observed change in freeze- 
out phase space density |^ ^ . Finally, we turn to the 
dependence of the scattering rate and escape probability 
on temperature. In doing so, we find that « 50% of the 
particles are emitted at temperatures below ~ 100 MeV. 
We discuss this result, together with other findings of our 
calculations in Section IVI 



II. THE ESCAPE PROBABILITY 

The following discussion of the freeze-out process 0, 
3 is based on the particle escape probability V{x,p,t) 
defined as 



/osc(a;,p,T) ^ r{x,p,T)f{x,p,T) . 



(1) 



Here, /(x,p, r) is the distribution function of a given par- 
ticle species and /esc(a;,p, t) denotes the distribution of 
that fraction of particles which have decoupled from the 
system prior to time t and will not rescatter in the fu- 
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ture. For a particle with momentum p at space-time 
point {x,t), the probability to escape the medium with- 
out future interaction is [3, Q 



T'{x,p, t) = exp 



drTZ^x + vT,p) 



(2) 



Here, the scattering rate TZ{x,p) denotes the inverse of 
the mean time between collisions for a particle at posi- 
tion X with momentum p. The opacity integral in ^ 
determines the average number of collisions of the parti- 
cle after time r. Thus, for a particle which tries to escape 
the medium at time t with velocity v, eq. ^ determines 
the probability that the particle succeeds to do so. 

The limiting case of a sharp freeze-out along the three- 
dimensional proper time hypersurface given by Tfi- Q can 
be characterized by an abrupt change of the particle es- 
cape probability (O from zero to one 



'P{x,p. 
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T < Tfi- 
T > Tfr 



(3) 



Freeze-out along more general three-dimensional hyper- 
surfaces can be defined by requiring this kind of 
threshold behavior of the escape probability along those 
hypersurfaces. It corresponds to a scattering rate which 
changes abruptly from a large value to zero along the 
freeze-out hypersurface. 

Hypersurfaces of the type Q are characterized by a 
criterion which tests only the medium, but does not de- 
pend on the particle momentum and particle properties. 
In the general case of continuous freeze-out, the situation 
is different since freeze-out does depend on particle mo- 
mentum and properties. To characterize the four- volume 
from which particles decouple, one requires Q 



p^d^V{x,p,T) > 0. 



(4) 



This specifies the region in which the amount of escaped 
particles of momentum p^^ is growing. Note that this gen- 
eral condition also characterizes the freeze-out hypersur- 
face in case of the sharp freeze-out according to eq. Q. 
Condition ^ assumes, however, that there is no addi- 
tional particle production, i.e., p^d^f{x,p,T) — 0. In 
the presence of such particle production, one generalizes 
(El) to d^ipi'Vf) > 0. 

In general, the evaluation of the criterion is com- 
plicated. In the Discussion, we shall consider a related, 
much simpler condition which gives some access to the 
structure of a four-dimensional freeze-out region on the 
basis of the escape probability. 



III. THE DEPENDENCE OF PARTICLE 
FREEZE-OUT ON EXPANSION 

It has been argued a long time ago Q that freeze-out 
happens when the dilution rate becomes larger than the 
collision rate. Here, we demonstrate how the collision 



rate and the expansion strength determine the escape 
probability. 

To this end, we consider the simplified case of a par- 
ticle of vanishing transverse momentum in the center of 
the fireball. Such a particle does not propagate through 
layers of different density but — due to expansion — finds 
itself in a medium of decreasing density. This simplifica- 
tion allows to illustrate the effect of dynamical properties 
of the collision region on freeze-out without being sensi- 
tive to further complications of the general case such as 
finite-size effects which depend on the production point 
and velocity of the test particle. We start the discus- 
sion of the escape probability (|2J) by considering a simple 
power-law ansatz for the scattering rate 



n{p = 0, r = 0, r) = 7^(T) = 7^o 



To 



a>l, (5) 



Here, TZq is the scattering rate at the time tq. Below, we 
detail soon the assumptions on which Q is based and 
how the exponent a characterizes the expansion strength 
of the system. To set the stage, however, we consider first 
the opacity integral of Q, 



/"OO 

/ dT7^(r) 

J Tn 



a ~ 1 



(6) 



This illustrates the typical interplay of scattering and 
expansion: a given escape probability can be obtained 
for different values of the scattering rate, with a higher 
scattering rate TZq compensated by a stronger expansion 
(larger a). 

We further illustrate this point with two analytically 
accessible models for the fireball expansion. Both are 
based on a factorized ansatz for the scattering rate in 
terms of the (averaged) cross-section a for scattering in 
the medium, and the average velocity Vre\ relative to 
other particles. 



TZ{t) = ap{T)vi.ci . 



(7) 



While the time dependence of TZ does not factorize in 
the general case, one may hope to capture the dominant 
features of a realistic dynamical evolution by retaining 
the time dependence of the density, only. 

For the first model, we choose a power-law fall-off of 
the density 



Pir) 



Po 



(8) 



where po is the density at the time tq. From this, the 
expression © can be recovered. This time-dependence 
of the density may not be realised during the whole evo- 
lution of the fireball; we assume it just in the final stage 
of the collision. Therefore, the time r is now defined by 
the dilution rate 7 



1 dp 
p dr 



a 

T 



(9) 
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Even though we are interested in the situation along 
the longitudinal symmetry axis of the fireball, the time 
r agrees with the "usual" longitudinal proper time 
\Jt'^ — only if the density evolves according to (|SJ) from 
the very beginning of the collision (i = 0). Otherwise it 
corresponds to a different starting point to of the time 
scale I3II 



(10) 



This is analogous to the Hubble time in cosmology which 
is defined as the inverse of the expansion velocity gradi- 
ent. We insert this dynamical information in the usual 
parametrization of the expansion four-velocity 

= (cosh 77 cosh 774 (r, r) , cos sinh r/t (r, r) , (11) 
sin0 sinh774(T, r) , sinh?] coshryt(T, r)) . 

In accord with relation (|10|l . the space-time rapidity is 
defined as 



1 



77 



■In 



{t -to) + z 



[{t-to)-z 



and r, (jj are the standard radial coordinates used in the 
plane transverse to the beam. The transverse rapidity 
r]t will be assumed to grow linearly with the radial coor- 
dinate. If the particle number is conserved, which is a 
good assumption at the end of the hadronic phase, the 
dilution rate is related to the divergence of the velocity 
field iia] 



Idp 
p dr 



1 

P 



= — u''d^p = d^u^. 



(12) 



If we assume that rjt{T,r) (x r, then eq. (|12|l links the 
four- velocity field Hll() with the ansatz ((SJ. This specifies 
the time dependence of the transverse rapidity 



T 

a — 1 



(13a) 
(13b) 



In the vicinity of r = 0, the choice x = 1 leads to quasi- 
inertial flow which corresponds to an asymptotic solution 
of the fireball hydrodynamics The values % < 1 

and X > 1 stand for radially decelerating and accelerat- 
ing flow profiles, respectively. This illustrates the phe- 
nomenological consequences of a specific choice of the 
exponent a. Relation (|13b|l allows to rewrite the opacity 
integral © 



flow and still lead to the same escape probability. We 
note that it is not the average flow velocity, but the local 
flow gradient which determines the local density decrease 
and thus enters (|14ll . In Appendix 1X1 we relate this result 
to the freeze-out criterion requiring the scattering rate to 
be smaller than the dilution rate |3] . 

For a second simple dynamical model of the fireball 
evolution, we turn to recent hydrodynamic simulations 
. These indicate that the transverse expansion at the 
freeze-out stage may be better described by the ansatz 



Vt{T,r) = £,r, 
^(t) = const , 



(15a) 
(15b) 



rather than by the expressions (jHl) and H13(l used above. 
In this case, the variable r is still given by the relation 
(|10|l . but corresponds now to the inverse longitudinal gra- 
dient of the expansion velocity 



3\-l 



The ansatz (jl5|l implies a different time-dependence of 
the density. From eq. H12|) . one finds for the density at 
the center of the fireball 



p{t) = Po— exp {-2£_{t - tq)) 



(16) 



where po is again the density at the time tq. The corre- 
sponding opacity integral reads 



/•OO 

/ dr 7^(T) = 7^o TO exp(2^ro) r(0, 2^to) , 

J To 

where r(a, x) is the incomplete F-function ^15.] 



(17) 



T{a,x) ~ 



dte-^e-^ . 



= r(a) - 



dtf-'^e-*. 



Relation H17|) is less transparent than H14I) but it shows 
the same qualitative feature: stronger scattering can be 
compensated by transverse expansion and lead to the 
same escape probability. 



IV. THE SCATTERING RATE IN A THERMAL 
MODEL 



dTTZiT) 



a — 1 



no 
2 



(14) 



According to (|13a|l . x/'Hd is the gradient of transverse 
rapidity at the time tq. At r = 0, it is to a good approx- 
imation equal to the transverse velocity gradient. This 
substantiates the statement made above: a higher scat- 
tering rate can be compensated by a stronger transverse 



In the previous section, we studied the effect of ex- 
pansion on freeze-out. We derived expressions for the 
opacity integral which depend on the scattering rate at a 
fixed time tq. Here we calculate this scattering rate and 
we study how it depends on hadrochemistry and tem- 
perature. In particular, we determine the scattering rate 
corresponding to the hadronic final states at SPS and 
RHIC. 
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The scattering rate for a test pion of momentum p due 
to interactions with particles of type i and momentum k 
can be written as [m. Il7l| 



2^m\ + p'^y^mf + fc2 



Here, cri(s) denotes the total cross-section for collinear 
collision, 



Sa = (mj + m^)^, 



(19) 
(20) 



and the cms energy s = s{k,p). 

We assume that the distribution of scattering partners 
Pi(x,k) is given by the equilibrium form. In this way 
we neglect the modification of the distribution function 
due to decoupling of some particles. Such approximation 
was argued to cause only a small error on the escape 
probability in the practically relevant cases Q ■ Thus we 
write 



Pi{x,k) 



9i 



(21) 



where Qi is the degeneracy of the species, T is the tem- 
perature, Ek = ■\/fc2 + raf and the chemical potential jjLi 
fixes the total density. This is known to provide a good 
description of the hadronic final state at SPS and RHIC 

Mm 

The total scattering rate is obtained by integrating the 
expression ()18|l over momentum k and summing over all 
species i 



n{x,p) = J2 J d'k 



dRi{x,p) 
'¥k 



(22) 



All pion scattering rates are computed in the rest frame 
of the hadron gas. We include pions, (anti)nucleons, 
kaons, rhos, and (anti)deltas as scattering partners. The 
total cross section of pion-baryon scatterings is parame- 
terized as ,20;1 



j^, m^\\Jr, Mr) {2Sr + 1) 
{2S, + l){2S^ + l) 



>Ls {Mr-Vs?+TlJ^ 



(23) 



This is the usual Breit-Wigner resonance formula where 
Mr is the resonance mass and Ftot and Tr—>-ni are the 
total and the partial width for the given decay chan- 
nel, respectively. Summation in 123|l runs over all rele- 
vant resonance states listed in (2fli] . Momentum in the 
cms is denoted as Pcms- The pre- factor takes care of 
proper counting of spin and isospin states. The square 
of Clebsch-Gordan coefficient (ji, rrii, j^, \ \Jr, Mr) 
assures that the appropriate fraction of the resonance 
state is picked in coupling of the isospin states of the two 
scattering partners. Particle properties are taken from 



|2l| except for higher excitations with large uncertain- 
ties which are taken from [23|- For scattering on both 
strange and non-strange mesons, the formula ^2',\\ is used 
again and a 5 mb momentum-independent contribution is 
added to account for elastic processes. The cross-section 
with baryons is taken to saturate at 30 mb for above 
3 GeV/c. This, however, has a negligible effect on the cal- 
culation since high momentum particles are suppressed 
in the thermal distribution. 

To fix the hadrochcmical composition of scattering 
partners at SPS and RHIC, an estimate of the chemical 
potentials entering (|21|l is needed. We take the chemical 
potential for direct pions from data on pion phase-space 
densities, and those for all other particle species from ra- 
tios of dN/dy at midrapidity. For SPS at = 17^GeV 
and RHIC at = 130 ^GeV we determine the chemical 
potentials for three different temperatures: 90, 100 and 
120 MeV. 

In more detail: the pion chemical potential at the SPS 
is obtained by comparing the m_L-dependence of the av- 
erage phase-space density measured by NA44 and 
WA98 23] to the result expected from a thermalized 
boost-invariant source with box transverse density pro- 
file and a transverse flow profile rjt = v^jy/r/iibox [24| . 
The transverse momentum spectra relate transverse flow 
and temperature, such that -qf ~ 0.7 corresponds to T = 
90 - 100 MeV and 77/ ~ 0.55 for T = 120 MeV. From 
the ratios of dN/dy at mid-rapidity measured for pi- 
ons, protons, antiprotons, positive and negative kaons 
r25; 26J, particle densities and corresponding chemical 
potentials are extracted under the assumption that all 
particles originate from the same thermal source, in ac- 
cord with the use of equilibrated distribution in eq. I|18f) . 
Resonance decay contribution to pion production is ac- 
counted for. The density of neutrons is assumed to be 
the same as that of protons. Chemical potentials for p's, 
A's and A's are deduced by requiring detailed balance 

m 

^Ip = 2^^ , HA = Pp + P-TT , fJ-A = fJ-p + Mtt ■ (24) 

The chemical potentials extracted in this way are sum- 
marized in Table The upper and lower values of /i's 
were chosen such that our parameters over- or under- 
predict the data by at least the maximal amount allowed 
by the quoted error bars. This ensures insensitivity of 
our conclusions against fine-tuning of parameters. 

The same procedure is repeated for RHIC at y/s — 
130 AGeV. Results are summarized in Table ^1 The 
yields at mid-rapidity for pions, (anti)protons, and kaons 
were taken from 28] . Unlike at SPS, the feed-down to 
(anti)proton yields from weak decays was not corrected 
for in the used data. We performed a simple correction by 
assuming that the yields of A's and S's are given by the 
same temperature and their chemical potentials equal to 
that of the protons. The preliminary data on pion phase- 
space density presented recently by STAR show small 
statistical error bars only and lie close to the u ppe r 
bound of the large systematic errors quoted earlier 
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temperature 


90 MeV 


100 MeV 


120 MeV 




477 - 497 


435 - 455 


350 - 379 


Hp 


301 - 321 


238 - 259 


114 - 143 


fJ-TT 


50 - 65 


38 - 53 


10 - 30 


MA 


527 - 562 


473 - 508 


360 - 409 


MA 


351 - 386 


276 - 312 


124 - 173 


Hp 


100 - 130 


76 - 106 


20 - 60 


fJ-K 


162 - 182 


132 - 153 


72 - 102 


Hk 


111 - 131 


76 - 96 


4-33 



TABLE I: Chemical potentials in units of MeV used in the 
calculation for SPS at ^/s = 17 AGeV. 



temperature 


90 MeV 


100 MeV 


120 MeV 


Hp 


442 - 489 


393 - 452 


291 - 373 


Hp 


407 - 458 


354 - 418 


245 - 332 


Htt 


78 - 100 


70 - 97 


50 - 85 


MA 


520 - 589 


463 - 549 


341 - 458 


Ma 


485 - 558 


424 - 515 


295 - 417 


Hp 


156 - 200 


140 - 194 


100 - 170 


HK 


194 - 242 


170 - 229 


120 - 202 


Hk 


179 - 230 


153 - 217 


100 - 187 



TABLE II: Chemical potentials in units of MeV used in the 
calculation for RHIC at ^ = 130 y4GeV. 



Our analysis accounts for this not fully clarified experi- 
mental situation byassociating to the central value of the 
most recent data Q the large (asymmetric) systematic 
errors of To avoid a bias in relating SPS and RHIC 
energies, we compare the upper bound of the chemical 
potential for RHIC in Tabled (which corresponds to the 
central value in 9]) to the upper bound listed in Table HI 
for SPS. We thus underestimate the difference between 
RHIC and SPS. 

Figure shows the corresponding scattering rates for 
T = 100 MeV. The smaller contribution of nucleons at 
RHIC is largely compensated by antinucleons. Despite 
the increased pion phase-space density at RHIC, the pi- 
onic contribution does not dominate the total scattering 
rate because of the small pion-pion cross section. Scatter- 
ing on pions is of comparable importance to scattering on 
nucleons and antinucleons: both particle species lead to 
a momentum averaged scattering rate of ~ 0.2 (fm/c)~^. 

We observe that the increase of pion phase-space den- 
sity at RHIC when compared to SPS l^jllfl '^o signifi- 
cant impact on the scattering rate. The earlier suggestion 
that the pion phase-space density at the freeze-out should 
be a universal quantity Q did not account for the con- 
tribution of other particle species to the pion scattering 
rate. 

Figure |21 summarizes our results for the various sets 
of freeze-out temperature and chemical potentials which 
limit the range of values consistent with data from SPS 
and RHIC. Although particle densities are approximately 
the same for all temperatures, the scattering rate in- 
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FIG. 1: Pion scattering rate as a function of pion momentum 
at T = 100 MeV and the highest possible values of chemical 
potentials for SPS (left column) and RHIC (right column). 
Contributions to the total scattering rate from scattering on 
nucleons, antinucleons and pions are indicated. The lower row 
shows the baryonic and mesonic relative contributions. 

creases significantly with temperature. On the other 
hand we see a generic decrease of the scattering rate with 
the pion momentum. This behavior is mainly dictated by 
phase-space available for the collisions, but also depends 
on the cross-section as a function of s. We discuss this 
behavior in more detail in Appendix IbI 

In the lowest panel of Fig |21 we compare the extreme 
results for SPS and RHIC. If the scattering rate in one 
of the systems was clearly larger than in the other one, 
this would indicate a stronger transverse expansion in 
that system, as follows from Section HTll We see that the 
present data do not allow to make any conclusion on this 
subject just by using the scattering rate. 

V. CONCLUSIONS 

In this paper, we have have calculated the scatter- 
ing rates of pions which are characteristic for the late 
hadronic phase of the collisions at SPS and RHIC. We 
also illustrated for simple examples how these scattering 
rates determine the escape probability, and how the es- 
cape probability is affected by velocity gradients. More- 
over, we found that the temperature dependence of the 
scattering rate allows to constrain the range of tempera- 
tures at which pion emission is significant. To estimate 
this temperature range, we combined the scattering rates 
of Section Hvl with the expressions for the opacity inte- 
gral (|14|1 and (|17|l . Here, we discuss what is needed to 
overcome the limitation of our approach and what our 
results imply for freeze-out in the realistic case. 

To go beyond the calculations presented here will re- 
quire a refined dynamical model which follows the full 
space-time evolution of the fireball. Most likely, this will 
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FIG. 2: Pion scattering rate as a function of pion momentum 
calculated for temperatures 90, 100, and 120 MeV and all 
corresponding sets of chemical potentials from Tables |l] (SPS 
at ^ = 17 AGeV) and|ITl(RHIC at = 130 ^GeV). Lowest 
panel summarizes the lowest and the highest curves from both 
SPS and RHIC. 



involve numerical simulations. In particular, one may 
hope to extend in this way the expressions for the opac- 
ity integral in Section IIIII to the case of particles with 
finite momentum. For a realistic scenario, we envisage 
the following additional effects: 

(i) In a transversely expanding fireball, particles with 
p± > are produced from parts of the system that 
co-move transversely. Their momentum relative to 
the medium p^^ is thus smaller than p± . The corre- 
sponding scattering rate TZlp^^) is thus higher than 
the naively expected value at p±. Hence, this effect 
lowers the escape probability. 

(ii) A particle with finite p± can escape the fireball re- 
gion in a finite time. Thus, the corresponding opac- 
ity integral |(2Jl receives a contribution from a finite 



time only, in contrast to the infinite time which is 
accumulated by a test particle with pj_ = in our 
simplified calculation. This effect leads to a lower 
value of the opacity integral and enhances the es- 
cape probability. 

Also, the evaluation of the general freeze-out criterion 
(@J) is complicated and requires a complete knowledge of 
the space-time evolution of the fireball. To gain some 
understanding of the freeze-out process from the mod- 
els studied here, however, it is instructive to turn to 
a simpler condition. Let us define freeze-out for parti- 
cles at position x with momentum p to occur if the cor- 
responding escape probability V{x,p) increases above a 
certain threshold value. Doing this, we actually specify a 
three-dimensional hypersurface on which a given fraction 
/esc// = V oi all particles is already decoupled. How- 
ever, in contrast to previous freeze-out criteria [7». i8i. illl| . 
this condition does not only depend on the medium, but 
also on the particle momentum. Moreover, we can ac- 
cess some properties of the whole four-dimensional freeze- 
out region by varying the threshold value for V. In 
this way, the above model calculations give some insight 
into freeze-out criteria which are more general than the 
Cooper-Frye type. To put in numbers: A realistic esti- 
mate of the transverse flow gradient ^ = 0.08 fm^^ and 
the Bjorken freeze-out time tq — lOfm/c is taken from 
[l8l |. The value ^ = 0.08 fm~^ is also in agreement with 
hydrodynamic simulations |T^ . In order to achieve the 
escape probability of at least 1/3, the opacity integral 
must be smaller than — In 1/3 ~ 1.1. Then, eq. I|14|l 
leads to a scattering rate of less than 0.18 (fm/c)^^ while 
eq. (|17() puts the upper limit for the scattering rate to 
0.26 (fm/c)~^. At a temperature of 120 MeV, such a low 
value of the scattering rate is possibly reached at RHIC 
for p„ > 0.4GeV/c and at SPS at p^ > 0.3 GeV/c, as 
seen in Fig. [2| If we increase the required escape proba- 
bility to 1/2, the scattering rate must be smaller 0.11 or 
0.16, according to eqs. H14|l and (|17|l . respectively. For 
particles with < 0.4 GeV/c, the temperature has to 
drop to 100 MeV if V should reach 1/2. Particles with 
Ptt < 0.2 GeV/c must wait even longer. On the other 
hand, if we put the required value of V to 0.1, eq. H17|) 
leads to an upper bound of 0.54 (fm/c)~^ for the scatter- 
ing rate. In the case of T = 120 MeV and high /x's at 
RHIC pions with > 0.25GeV/c fulfill this condition; 
in the less extreme cases, the 10% escape probability is 
reached by almost all pions. 

This illustrates that particles decouple from the fireball 
gradually. It indicates that w 10% of the pions are de- 
coupled at T = 120 MeV, but approximately half of them 
will escape at local temperatures below 100 MeV. Thus, 
a large fraction of the particles decouples at rather low 
temperatures. We note that a temperature of 120 MeV 
was assumed as a freeze-out temperature in a recent work 
where a pion mean free path « 1 fm was suggested as 
the universal freeze-out criterion. This work neglects the 
temperature dependence of the escape probability and 
its dependence on the momentum of the particle. Also, 
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it does not consider the possible cancellation between 
stronger scattering and stronger transverse expansion. In 
our approach, a thermal distribution of T = 120 MeV at 
RHIC results in a pion mean free path of only 1.7 fm, 
only slightly larger than JJ. However, only 10 % of all 
particles are decoupled at this temperature. At lower 
temperatures, where the fraction of decoupled particles 
reaches 50%, we find a significantly larger average pion 
mean free path of at least 3-5 fm. 

Let us further comment on the significance of the mo- 
mentum dependence of the scattering rate. As discussed 
earlier, at certain p± the particle has a momentum with 
respect to the surrounding medium < P±, and thus 
the p_L-dependence of the scattering rate can be flatter 
than what is plotted in Fig. 13 We expect, however, that 
the monotonic increase of the escape probability with 
p± will be robust against refinements of our calculation. 
This suggests that in contrast to the idealization of a 
sharp freeze-out along a three-dimensional hypersurface, 
particle escape is ordered in momentum with high-px 
particles freezing-out earlier and thus originating from 
a smaller fireball. This effect could lead to a stronger 
transverse mass dependence of HBT radii 5^ than what 
is found in current model studies based on a Cooper-Frye 
freeze-out condition. These latter models (see e.g. the 
blast-wave model discussed in "so]) have difficulties in re- 
producing the observed strong Afx-dependencc of HBT 
radii at RHIC. In what concerns the mechanism proposed 
here, one may envisage compensating effects: for exam- 
ple, if the transverse ffow gradient is stronger at later 
times, then this will reduce the homogeneity regions mea- 
sured by HBT radii at later times. Such a scenario seems 
unlikely since it implies acceleration of the collective ex- 
pansion at a time when particles slowly cease to inter- 
act. However, to substantiate this expectation, a model 
which fully implements the fireball dynamics is needed. 
In our opinion, this warrants further investigations. In 
particular, hydrodynamical simulations should be revis- 
ited, which are typically based on sharp 3-dimensional 
freeze-out and thus miss any contribution from freeze- 
out along finite four- volumes. 

We thank J.G. Cramer, U. Heinz, J. Pisiit, and 
K. Redlich for fruitful discussions. 



APPENDIX A: INTERPLAY OF THE 
SCATTERING RATE AND THE DILUTION 
RATE 

Here, we give further details on the relation between 
the dilution rate and the scattering rate. We start from 
the opacity integral lO and we express the exponent a 
in terms of the dilution rate 7 via eq. © . This leads to 



/•oo 

/ ^^r7^(r) = 

J To 



a-1 



7- 



Tlo_ 
j_ 

To 



(Al) 



This can be directly compared to the statement that 
freeze-out happens when the dilution rate is at least as 



big as the scattering rate. For the scenario of continuous 
freeze-out considered here, a reasonable escape probabil- 
ity is reached when the opacity integral is of order one 
or smaller. In the model of eqs. ||SJ) and ([7|) this leads to 
the condition 



To 



i.e., there is a factor — I/tq in addition to the standard 
criterion 7 > TZq Q. The details of the relation between 
the scattering rate and the dilution rate at the freeze-out 
will depend on the particular time evolution of TZ and 
p, but the feature that larger TZ can be compensated by 
larger 7 is generic. 



APPENDIX B: GENERAL EXPRESSION FOR 
THE SCATTERING RATE 

In this appendix we derive a general expression for the 
scattering rate. This allows us to investigate its depen- 
dences on temperature and the test particle momentum. 

Following eq. 1)22(1 we write the scattering for some 
particle species of mass m. Assuming a Boltzmann dis- 
tribution of the scattering partners and suppressing the 
position dependence. 



nip) 



(Pk 



■g exp 



X (t(s) 



T 



\/{s- Sa){s - Sb) 

2 + p'^\^m'^ + k 



(Bl) 



If we change the integration variables to the center of 
mass energy s and the energy of the scattering partner 
Ek, and perform the i?fe-integration, this integral leads 
to 



7^(p) 



Sir'^py/p'^ + 



X / dsa{s)\J {s — — m?Y ^ ^ni^m? 

( P 



X sinh 



V2rm2 



Vi 



2 — T-vi 2 \ 2 



s — rat — m 



4™?.? 



xexp(-^^^±3(s-mj-m^)) . (B2) 



This scattering rate can be expressed as a convolution 
of the cross-section with the "distribution of two-particle 
approaches" "Dip, s) 



n{p) = / V{p, s) a{s) ds . (B3) 



8 



By comparing with (|B2|) . we obtain 



V{p, s) 



X sinh (^ 2j^ 2 ^^^ ^ ~ w?Y ~ 4rTi.^m" 



X exp 



In the arguments of the exponential functions of (|B4p . 
the momentum of the test particle and the temper- 
ature are inversely related. Thus if one ignores the 
l/p\/p^ + ra^ pre-factor, increasing T has the same ef- 
fect as decreasing p. This offers an explanation why 
TZ increases with raising the temperature but decreases 
when the momentum becomes larger. There is no colli- 



sion dynamics going into the calculation of 2?; it repre- 
sents merely the phase-space populated with thermally 
distributed scattering partners. 

However, the temperature and momentum dependence 
of the scattering rate is not given solely by the distribu- 
tion of two-particle approaches, but also reflects details 
of the convolution of s) with the scattering cross sec- 
tion cr(s). This can be seen e.g. by comparing our results 
with those obtained in [T^ . In that paper, the authors 
used a simple prescription a{s) ca 5{s — Af^) which leads 
to a qualitatively different dependence of the scattering 
rate on the test particle momentum. We verified that 
this is a consequence of neglecting the width of the reso- 
nance. By only assuming the lowest resonance states in 
a given channel and narrowing the resonance shape, we 
were able to reproduce the results of |l6j | as a limiting 
case of our calculation. 
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Technically, if we were solving for the whole fireball evolu- 
tion, we would have to match the prescription (|SJ| and its 
first time-derivative to the time-dependence of the den- 
sity at earlier times. This is done by tuning po and tq. 
Therefore, tq has to be a free parameter, and is not nec- 
essarily equal to the longitudinal proper time measured 
from the moment of the first approach of the colliding 
nuclei. 



